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Extinction law in ultraluminous infrared galaxies at z~l 
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ABSTRACT 

We analyze the multi- wavelength photometric and spectroscopic data of 12 ultra- 
luminous infrared galaxies (ULIRGs) at z '^ 1 and compare them with models of stars 
and dust in order to study the extinction law and star formation in young infrared 
(IR) galaxies. Five extinction curves, namely, the Milky Way (MW), the pseudo MW 
which is MW-like without the 2175 A feature, the Calzetti, and two SN dust curves, 
are applied to the data, by combining with various dust distributions, namely, the 
uniform dust screen, the clumpy dust screen, the internal dust geometry, and the 
composite geometry with a combination of dust screen and internal dust. Employing 
a minimum x^ method, we find that the foreground dust screen geometry, especially 
combined with the 8-40 M(^ SN extinction curve, provides a good approximation 
to the real dust geometry, whereas internal dust is only significant in 2 galaxies. The 
SN extinction curves, which are fiatter than the others, reproduce the data of 8(67%) 
galaxies better. Dust masses are estimated to be in excess of ~ IO^Mq. Inferred ages 
of the galaxies are very young, 8 of which range from 10 to 650 Myr. The SN-origin 
dust is the most plausible to account for the vast amount of dust masses and the fiat 
slope of the observed extinction law. The inferred dust mass per SN ranges from 0.01 
to 0.4 Mq/SN. 
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1 INTRODUCTION 

Recent developments in far- infrared (IR) and submillimeter 
astronomy have revealed new populations of galaxies associ- 
ated with thermal dust emission, such as dusty quasars and 
infrared galaxies. Thermal dust emission detected in quasars 
and galaxies at high redshift exhibits inferred dust mass in 
excess of ^ 10^M(^. Such a vast amount of dust observed at 
z > 6 (|Bertoldi et al .11 20031 : [Robson et al.]l2004llBeelen et al.1 
I2OO6I I strongly suggests that the dust has been condensed 
in ejecta of core- collapsed supernovae (SNe II) promptly af- 
ter the onset of the initial star formation. Asymptotic giant 
branch (AGB) stars are known to be the major source of 
dust in the present-day universe. However, the importance 
of AGB stars at high redshift is in controversy, because they 
inject dust grains into the ISM only after their progenitor 
stars have evolved off the main sequence. While a delay time 
of > 500 M y r since their formation has been suggested (e.g., 
lDweklll998l \ IValiante et al.1 (|2009l ) proposed that AGB stars 
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can dominate dust production in a delay time of > 1 5 Myr 
if the starburst is instantaneous. iMichalowski et al.l (|201Q| ) 
pointed out that AGB stars are not efficient enough to form 
dust in dusty quasars at z > 5, whereas SNe II may be able 
to account for the dust, and underlined the importance of 
dust grain growth in the interstellar medium (ISM). 

Quasars have been frequently used to study the ex- 
tinction law at high redshift. This is because of their sim- 
ple dust geometry, in which the central point-like emitting 
source is surrounded by dust grains which are physically sep- 
arated from the source itself. This is analogous to the clas- 
sical configuration assumed when dereddening single stars. 
There is a line of evidence that active galactic nuclei (AGNs) 
have extinction curves that differ from those observed in 
the Milky Way ( MW), Large and Sm all Magellanic Clouds 
(LMC and SMCl. lGaskell et al.1 (l2004l ) analyzed a large sam- 
ple of AGNs and found indications of a very flat extinction 
curve in the ultra-violet (UV). On the other hand, exploring 
large samp les of quasars f r om t he Sl oan Digital Sky Sur - 
vevfSDSS) jRichards et all (|2003l ) and lHopkins et al.1 (|2004l ) 
found that the extinction law in quasars at z < 2.2 is de- 
scribed by SMC reddening but neither by LMC and SMC 
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reddening nor by that found in the Gaskeh et alJ ((200^) sam- 
ple. Recently, iGallerani et alJ (|2010|) proposed that quasars 
at 2 = 3.9—6.4 have reddening that deviates fro m that of the 
SMC with a tendency flatten at A < 2000 A. Maiolino et al.1 
(l2004l ) found that SDSS J1048+4647, a broad absorption 
line quasar (BALQSO) at z—Q.2^ requires a flat extinction 
curve which is quite different from those observed in the 
MW, LMC and SM C, but in excellen t agjree ment with SN 
extinction curves bv lTodini fc Ferraral (|200lh . 

Although it is reasonable to expect that the extinction 
law in young galaxies would be dominated by SN-condensed 
dust, the complications introduced by dust geometries have 
hampered this line of studies; the attenuation of extended 
starlight by dust depends not only on the optical properties 
of dust b ut als o on the distribution geometry of the dust. 
ICalzettil (|200lh derived the empirical extinction law in lo- 
cal starbursts and blue compact galaxies, assuming the dust 
geometry described in a foreground screen, i.e., an appar- 
ent attenuation of starlight to be e"^"" . If a dust geometry 
in their galaxies is indeed a foreground screen, the slope 
(i.e., wavelength dependence) of the extinction curve (i.e., 
the optical depth Tt^) is significantly flatter than those ob- 
served in the MW, LMC, and SMC. The SED and spec- 
troscopic data of a young, ultralum inous infrared galax y 
(ULIRGJj dX z r^ 1 were analysed bv lKawara et al.l (|201lh , 
who found that the data require a flat extinction curve and 
a dust geometry of a foreground dust screen, resulting in 
better agreement with the extinction curve sy nthesized from 
SN-condensed dust bv lHirashita et al.l (|2005l ) than those ob- 
served in local starbursts, MW, LMC, and SMC. 

The slope of the extinction curve is expected to reflect 
the sites of dust formation and the subsequent processing in 
the ISM. The extinction curves can be ordered as the SMC, 
LMC, MW, and Calzetti curves from the steepest to the 
flattest. The synthetic extinction curves of newly condensed 
dust grains in SN ejecta before the ISM reprocessing (here- 
after called original SN extinction curves) are flatter than 
the Calzetti curve ( Todini fc Ferraral 200ll : iNozawa et ahl 



Table 1. Sample galaxies at z ~ 1 



I2OO3I : iMaiolino et all 1 20041: Iffirashita et"alT l2005l : see also 
Figure 3 in iKawara et al.l |201l|). Then, the bulk of 
dust grains condensed in the ejecta is destroyed through 
the passage of th e reverse shock (|Bianchi k, Schneiderl 



tne passage 01 tn e reverse snock (piancni <y bcnneiaen 
l2007llNozawa et al.1 l2007llHirashita et al.1 I2OO8I ) . Finallv. 



grains that survive the reverses shock are s hattered in 
the warm ionized ISM (Hirashit a et al.l |2010|). Nonethe- 
less, the final SN extinction curves, after taking into ac- 
count the reverse shock and shattering, are expected to 
have a slope as flat as those of the original extinction 
curves. The shattering may also be important in old galax- 
ies like MW, LMC, and SMC, where AGB stars are dom- 
inant sources of dust production, because some observa- 
tions (lGroenewegenlll997llGauger et al.lll999l ) suggest AGB 
stars form large grains. Dust grain growth in the ISM 
would be efficient in infrared galaxies where there must be 
many molecular clouds near SN remnants, and it is sug- 



^ Infrared galaxies are classified according to bolometric infrared 
luminosity Liji(8 — l,000|xm); Ljji > 10^^ Lq as luminous in- 
frared galaxies (LIRGs), Lju > 10^^ Lq as ultraluminous infrared 
galaxies (ULIRGs),a nd Ljr > 10^^ Lc;) as hype rluminous infrared 
galaxies (HyLIRGs) dSanders fc Mirabe]|[l99^ V 



No. 



Name 



Redshift Reference 



9 

10 
11 
12 



COSMOS J100121. 29+023535.8 0.850 

COSMOS J100125.12+024506.5 0.980 

COSMOS J100158.31+023505.2 1.021 

COSMOS J100143.83+015325.0 1.066 

COSMOS J100021. 79+014916.1 0.995 

COSMOS J100030.12+020513.3 0.894 

COSMOS J100123.22+021013.3 0.850 

COSMOS JlOOlOO.22+021503.6 0.893 

COSMOS J100008.36+022244.6 1.093 

COSMOS J095910.21+022523.3 1.026 

SST J1604+4304 1.135 

HRIO (ERO J164502+4626.4) 1.44 



2 

3,4 



Reference. - (1) Cosmos archive at http://irsa.ipac.caltech.edu: 
COSMOS Photometry Catalog January 2006: 
S-COSMOS IRAC 4-channel Photometry Catalog June 2007; 
S-COSMOS MIPS 24/im Photometry Catalog October 2008; 
S-COSMOS MIPS 70/im Photometry Catalog v3 January 2009; 
COSMOS MAGELLAN Data Version 3.1 (spectroscopy); 
zC OSMOS data release d i n 2008 October (sp ect roscopy); 



(2) iKawara et aIll201Ql : (3) IStern et al.ll2006l : (4) IPev et al.l [19991 . 

gested that the dust mass in the high- redshift quasars is 
dominated by grain growth in the IMS (|Dweklll998llDrainel 



l2009llMichalowski et al.ll2010llValiante et al.ll201lh . It seems 
that this process result in flat extinction curves in the sense 
that the larger grains show the flatter extinction curve. How- 
ever, the extinction laws found in the high-redshift quasars 
and local starbursts are flatter than those found in the MW 
for which IPraind (|2009|) attributes the bulk of the dust mass 
to the growth of grains in the ISM. 

As discussed so far, z < 2.2 quasars show a striking 
contrast in the slopes of the extinction curves against the 
z r^ 1 ULIRG and higher redshift quasars. In this paper, we 
explore this iss ue by expanding our previous work on a single 
young ULIRG (|Kawara et ahlboilh to 12 young ULIRGs at 
zr^l. 



2 DATA 

The sample consists of 12 infrared galaxies at 2; = 0.85—1.44. 
The redshifts and data sources of our 12 galaxies are 
shown in Table [D They are SST J1604+4304, HR10(ERO 
J164502+4626.4), and 10 galaxies found in the Cosmolog- 
ical Evolution Survey (COSMOS) field. We have searched 
archives and the literature for those galaxies which meet 
the following criteria: (1) the far-IR emission is detected; 
(2) the optical spectrum indicates a redshift of ^ 1, and 
gaseous emission lines such as [OII]A3727 and/or stellar ab- 
sorption lines such as Ca H/K and HI absorption lines were 
detected; (3) broad-band photometry has been done from 
the UV to mid-infrared; (4) there is no contamination from 
foreground galaxies; (5) the bolometric luminosity is dom- 
inated by the undergoing star formation with the criterion 
of /^(8|J,m)//^(4.5|J,m) < 1.9 (IPope et al.1 12008) and no ob- 
vious evidence of emission caused by an AGN such as X-ray 
emission. As a result, we have found these 12 galaxies all of 
which are classified as starburst-dominated ULIRGs. 
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Figure 1. Postage stamp images of the 10 COSMOS galaxies. The sizes of each images are 4^^x4^^ for F 814W to Ks, lO'^^xlO ^^ for 
3.6M-mto 8.0M,m, 2 0^^x20^^ for 24M,m, and 50^^x50^^ for TOj^m. North is up and east is to the left. Refer to lKawara et al.l (|2010l ) and 
IStern et al.l (|2006h for the images of No.ll (SST J1604+4304) and No.l2 (HRIO), respectively. 



Table 2. Optical and 2|j.m photometry of the 10 COSMOS galaxies 



No. 


u* 


g+ 


r+ 


i+ 


z+ 


Ks 


1 


0.530 ± 0.033 


0.924 ± 0.041 


3.04 ± 0.07 


8.43 ± 0.10 


14.2 ± 0.2 


62.6 ± 3.0 


2 


0.888 ± 0.048 


1.27 ± 0.05 


3.05 ± 0.07 


7.76 ± 0.11 


13.5 ± 0.2 


63.8 ± 2.9 


3 


0.096 ± 0.022 


0.164 ± 0.024 


0.656 ± 0.034 


1.95 ± 0.06 


3.78 ± 0.11 


38.0 ± 3.0 


4 


0.892 ± 0.052 


1.50 ± 0.07 


3.97 ± 0.09 


9.52 ± 0.14 


19.2 ± 0.2 


97.8 ± 4.6 


5 


2.20 ± 0.06 


2.51 ± 0.06 


5.10 ± 0.08 


9.46 ± 0.11 


13.8 ± 0.2 


44.0 ± 3.4 


6 


0.225 ± 0.023 


0.403 ± 0.003 


1.38 ± 0.05 


4.43 ± 0.08 


7.91 ± 0.13 


47.6 ± 2.9 


7 


2.07 ± 0.06 


3.22 ± 0.07 


6.50 ± 0.09 


12.6 ± 0.1 


18.4 ± 0.2 


49.3 ± 3.1 


8 


1.23 ± 0.05 


1.97 ± 0.06 


4.60 ± 0.09 


11.1 ± 0.1 


18.1 ± 0.2 


75.7 ± 3.5 


9 


1.10 ± 0.04 


1.26 ± 0.04 


2.13 ± 0.06 


4.03 ± 0.07 


6.98 ± 0.12 


30.8 ± 3.1 


10 


1.13 ± 0.04 


1.83 ± 0.05 


3.68 ± 0.07 


7.73 ± 0.10 


13.8 ± 0.2 


50.6 ± 3.0 



Note. - Total fluxes are given in units of /iJy. 
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Table 3. S'piteer photometry of the 10 COSMOS galaxies 



No. 


3.6M.m 


4.5|J.m 


5.8|J.m 


8.0M.m 


24M,m 


70|j,m 


1 


106 ± 0.2 


81.6 ± 0.3 


74.7 ± 1.1 


76.5 ± 2.8 


940 ± 252 


20900 ± 3300 


2 


89.6 ± 0.2 


73.1 ± 0.3 


66.0 ± 1.2 


58.9 ± 2.9 


797 ± 19 


11700 ± 2600 


3 


81.7 ± 0.2 


66.5 ± 0.3 


72.9 ± 1.2 


82.1 ± 3.0 


635 ± 20 


11100 ± 2600 


4 


128 ± 0.7 


111 ± 0.7 


78.6 ± 1.3 


87.9 ± 2.9 


780 ± 22 


11600 ± 2700 


5 


68.4 ± 0.2 


54.5 ± 0.3 


49.1 ± 1.2 


60.7 ± 2.8 


1000 ± 92 


14200 ± 2800 


6 


79.0 ± 0.2 


62.1 ± 0.3 


50.7 ± 1.2 


62.6 ± 2.7 


553 ± 16 


14900 ± 2700 


7 


67.2 ± 0.2 


50.1 ± 0.3 


62.5 ± 1.1 


81.4 ± 2.8 


1190 ± 256 


20800 ± 3100 


8 


93.0 ± 0.2 


69.8 ± 0.3 


68.9 ± 1.1 


56.4 ± 2.5 


849 ± 165 


9900 ± 2500 


9 


57.2 ± 0.5 


51.6 ± 0.5 


35.1 ± 1.2 


53.3 ± 2.7 


573 ± 179 


16100 ± 2300 


10 


91.2 ± 0.2 


71.3 ± 0.3 


46.5 ± 1.1 


53.4 ± 2.9 


590 ± 199 


7300 ± 1800 



Note. - Total fluxes are given in units of /xJy. 



The dat a for SST J1604+4304 were taken from 
iKawara et alJ (|2010h where the multi-wavelength data cov- 
ering the X-ray, UV/optical, infrared, and radio and the 
optical spectrum are given. For HRIO, IStern et alJ (|2006|) 



compiled the multi- wavelength data from the UV/optical to 
radio and the optical spectrum ( strictly speaking from 2750 
to 3850 A) is given in lPev et all ([^99). 

The photometric and spectroscopic data for the 10 
COSMOS galaxies are available in the COSMOS archive 
at [http://irsa.ipac.caltech.edu/data/COSMO S/ Figure [T]il- 
lustrates postage stamp images of the COSMOS galaxies 
from u* to 70|J.m. For full details of the ob servations of 
the C OSMOS galaxies, the reader refers to ICapak et al.l 
(|2007|) for optical to near-infrared imaging performed on 
the Hubble Space Telescope and ground based telescopes in- 
cluding the Subaru 8.3m telescope, the Kitt Peak National 
Observatory (KPNO) 4m telescope, the Cerro Tololo Inter- 
American (CTIO) 4m telescope, and the Canada-Fra nce- 
Hawaii (CFH) 3.6m telescope, to ISanders et all (l2007l) for 
the C OSMOS Spitzer survey (S- COSMOS).lLe Floc'h et al.1 
(I2OO9I ) for the 24|am observati ons, iFraver et al.1 (l2009l ) for the 



70 and 160 ^im observations, iTru mp et al](i20Q9l ) for 5600 



9200A spectrosco py on the the Magellan 6.5m telescope, and 
iLillv et al.1 (|2007l ) for 5500-9700A spectroscopy on the Eu- 
ropean Southern Observatory (ESO) 8m Very Large Tele- 
scope (VLT). 

For COSMOS galaxies, the optical counterparts of far- 
IR sources were selected as follows. We first selected objects 
at2; = 0.8 — 1.21in the spectroscopy catalogs and searched 
for the counterparts in the optical imaging data catalogs. If 
there is a single object in the optical imaging data within 
an error circle radius of 0.2^^ from the spectroscopic posi- 
tion, we identified the spectroscopic object with the object 
in the optical imaging catalogs. We then searched for the 
mid-infrared counterpart in the IRAC (3.6, 4.5, 5.6, S.O^tm) 
catalog. Again, if there is one IRAC object within an error 
circle radius of O.S^^from the position of the optical, spec- 
troscopic object, we regard the IRAC object as optically 
identified. This process was repeated for the 24|j,m counter- 
part within an error radius of 1.0^^ and 70|J,m counterpart 
within 4.0^^ We found no 160|J,m counterparts within 8.0^^ 
As a result, we found 10 star burst-d ominated ULIRGS at z 
= 0.8 - 1.2. lKartaltepe et al.l (|2010l l, applying more relaxed 
criteria than ours, found 50 70|am sources having spectro- 
scopic data in the same field and the same redshift range. 



Assuming that 70|am sources with > 10 mJy are ULIRGs at 
z > 0.8, 25 of the 50 sources are ULIRGs with spectroscopic 
data. They suggest that an AGN fraction is about 50%, thus 
starburst-dominated ULIRGs are about 12, fairly consistent 
with the number of ULIRGs in our sample. 

Optical data were taken in many photometric bands 
for COSMOS galaxies. If we use all these data in the fol- 
lowing analysis, the weight of the spectral energy distribu- 
tion (SED) in the optical will overwhelm the other data. 
We thus use u* taken on the CFH telescope, ^^,r^,'i^, z~^ 
on the Subaru telescope, Ks on the CTIO 4m telescope. 
Aperture photometry with a 3^Miameter is given in the cat- 
alog. We appl ied the corre c tion t o the total magnitude as 
instructed by ICapak et al.l (|2007h . The aperture photom- 
etry with a 2.9''diameter given in the S-COSMOS IRAC 
4-channel Catalog was converted to the total fluxes by ap- 
plying correction factors of 0.90, 0.90, 0.84, & 0.73 for 3.6, 
4.5, 5.8, & S.O^im, respectively (see README attached to 
this catalog). The total fluxes of the COSMOS galaxies are 
summarized in Tables [2] and [3] in units of /iJy. The SEDs of 
the 12 galaxies are plotted as a function of rest wavelength 
in Figure (2] 

The optical spectra of our 12 galaxies are illustrated in 
Figure [3] along with the vertical lines showing the positions 
of the [OII]A3727 emission line from HII regions and the 
stellar absorption lines, namely, MgIIA2798, Call H & K, H^, 
and He. [OII]A3727 indicates the presence of OB stars, while 
the absorption lines are useful t o measure the age of the 
relatively old stellar populations ([Kawara et al.ll2Qllh . The 
equivalent widths of the lines are given in Tabled along with 
the D(Call) index. This index is defined as D(Call) = 2 x 
EW(CaII K)/[EW(Hc)+EW(CaII H + H^)] (JKawara et al.1 
2010), where EW(CaII K), EW(Hc), and EW(CaII H + He) 
are the equivalent widths of Call K, H^, and the blend of 
Call H and He, respectively. 

Table [5] gives the 8-l,000^im luminosity Ljr and the 
mass of dust Md of the 12 sample ga laxies. We obtained 
LiR by fitting the lPale fc Heloul (120021 ) ULIRG template to 
the 70 |j,m data except for SST J1604+4304 for which the 70 
and 160 \im data were used for fitting. Deriving far-IR col- 
ors /i/(160^Lm)//iy(70|J,m) of 70|J,m sources that are not de- 
tected at 160M,m f r om th e stacking of the 160|am data points, 
iKartaltepe et al.l (|2010l ) found that those not detected at 
160|am have hotter far-IR colors than those detected at 160 
|j,m, and that the far-IR colors of 70 ^im sources not detected 
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Figure 2. 0.1 - 100 |Lim SEDs of the 12 galaxies, illustrating ufT^(ergs cm ^ s ■'^) as a function of rest wavelength(|j.ni). 



Table 4. Spectroscopic proterties of the 12 galaxies 



(1) 

No. 


(2) 
EW( Oil]) 

(A) 


(3) 
D(Call) 


(4) 
EW(Hc) 

(A) 


(5) 
EW(Ca K) 

(A) 


(6) 
EW(Ca H + He) 

(A) 


1 


45.6 ± 2.4 


0.28 ± 0.13 


10.6 ± 1.6 


3.5 ± 1.5 


14.2 ± 2.2 


2 


27.3 ± 2.7 


0.11 ± 0.13 


4.3 ± 1.9 


1.0 ± 1.1 


12.8 ± 1.5 


3 


10.1 ± 3.5 


0.27 ± 0.36 


12.6 ± 5.2 


2.6 ± 3.4 


6.8 ± 3.8 


4 


31.4 ± 3.3 


0.15 ± 0.05 


6.5 ± 1.4 


1.6 ± 0.6 


15.4 ± 1.2 


5 


65.3 ± 1.8 


0.18 ± 0.22 


6.4 ± 2.6 


1.0 ± 1.2 


4.7 ± 1.9 


6 


9.8 ± 1.8 


0.51 ± 0.26 


8.4 ± 3.7 


4.9 ± 2.2 


10.8 ± 2.6 


7 


35.0 ± 0.6 


0.13 ± 0.26 


1.8 ± 0.6 


0.3 ± 0.6 


2.9 ± 0.7 


8 


23.4 ± 1.3 


0.18 ± 0.41 


3.6 ± 1.0 


0.6 ± 1.2 


2.6 ± 2.3 


9 


63.2 ± 6.0 


0.02 ± 0.14 


6.1 ± 1.9 


0.1 ± 1.0 


9.0 ± 1.9 


10 


9.1 ± 2.1 


0.37 ± 0.26 


5.8 ± 2.0 


2.6 ± 1.8 


8.4 ± 2.4 


11« 


31.1 ± 4.3 


0.03 ± 0.11 


7.6 ± 2.1 


0.4 ± 1.4 


16.7 ± 2.2 


12^ 


47 ± 5 




EW(Mgll) = 7.6 ± 3.1 





Note. - Col 2 for the eqivalent width of the [ Oil] line at 3727A. Col 3 for the Call index defined as D(Call) =2 x 
EW(CaII K)/[EW(H^)+EW(CaII H + He)] (iKawara et al.ll2010l V Cols 4-5 for the equivalent width for the H^ fine 
at 3889A and the Call K line at 3934A, respectively. Col 6 for the equivalent width of the blends of the Call H & 

He lines at 3968 and 3970A, respectively.^ 

" The equivalent widths wer e taken fromlKawara et al.l ([20101) • 

b EW([OII ] ) was taken from lDev et al.l (|l999l Vand EW(MgIIA2798) was measured using the spectrum given by 
iDev et al.l (ll999l V 
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Figure 3. Optical spectra of the 12 galaxies. The relative flux densities (per A) are plotted as a function of rest-frame wavelength (A). 
The vertical lines mark the positions of Mgll at 2798A, [Oil] at 3727A, Call H & K at 3968 and 3934A respectively, H^ at 3889A, and H^ 
at 3970 A. 



at 160 M^m scatter around the IPale fc Heloul (|2002l ) ULIRG 
template at 2; '-^ 1. It is thus appropriate to utilize this tem- 
plate to derive Ljr for those not detected at 160 ^im. The 
typical uncertaint y for Ljr using; the 70|Li m data is estimated 
to be - 0.2 dex JKartaltepe et al.ll2010l ). SST J1604+4304 
was re- analyzed in the same way as the other galaxies. While 
the Arp 220 template was used to esti mate Lir and the dus t 
temperature in our previou s analy sis (jKawara et al.ll201ll ). 
we apply the IPale &: Heloul ^qM) ULIRG template in the 
present analysis. In addition, we include a new term of the 
dust mass in the minimum x^ method. The new results are 
quite consistent with our previous results. 

The dust mass was estimated using the formula of 
Md/U = Dl/[{1 + ^)A^^,B^,(Td)], where /^ is the observed 
flux density at frequency zy, Ue = {l-\-z)iy, Dl is the luminos- 
ity distance to the emitter, Bu^ (Td) is the Planck function 
for Td the temperature of dust at z^e, i^ue is the absorption 
cross-section per unit dust mass at i^e • As can be seen in this 
formula, Md depends on the product K,uBu{Td) oc v^ BuiTd) 
if Ku (X v^ . 

It is k nown that there i s an ant i- correlatio n between 
P and Td (JBlain et al.ll2003l : iBracco et al.ll201lh . This de- 
generacy may be due to the intrinsic properties of dust 
grains or application of a single-temperature model to the 



SEDs w hich average variou s temperatures along the line 
of sight (jBracco et al.ll201ll ). Thus, when fitting a single- 
temperature model to the SED, the resultant Td depends 
on P or the properties of the dust extinction curves; the ex- 
tinction curves we study here, are the MW, pseudo MW, 
Calzetti, SN curves, as fully described in section 3.1.2. 



The dust temperatures cannot be directly derived from 
our data, because of the lack of data lon^ward of 160 i^ni . 
We thus estimated Td by fitting to the lPale k Heloul (|2002l ) 
ULIRG template, resulting in Td = 40.5 ± 9.3 for the SN 
dust with P c^ 1.6 and Td = 37.0 ± 8.9 for the MW dust 
with P c^ 2] Td derived with the SN dust is 10% higher 
than that with the MW dust. This temperature difference 
is expected from the degeneracy between (5 and Td. Table 
E] gives Md for the 8-4OM0 SN dust. Md in Table [5] should 
be multiplied by 0.81 for the 8-3OM0 SN dust and 0.78 for 
the MW dust. It is noted that Md is derived by using n^^^ 
at 80 ^tm which is 44cm^ g~^ for MW dust,24cm^ g~^ for 
8-3OM0 S N dust, and 20cm^ r ~^ for 8-4OM0 SN dust (see 
Figure 7 in lKawara et al.ll201lh . 
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Table 5. Infrared luminosity, dust mass, and dust shell radius of 
the 12 galaxies 



No. 


LiR- (IQi^Lo) 


Md^ (lO^Mo) 


r^ (arcsec) 


1 


2.60 ± 1.14 


1.66 ± 1.06 


0.63 


2 


2.87 ± 1.25 


1.55 ± 1.02 


0.74 


3 


2.66 ± 1.18 


1.73 ± 1.15 


0.47 


4 


2.65 ± 1.20 


2.15 ± 1.43 


0.52 


5 


3.73 ± 1.63 


2.00 ± 1.30 


0.80 


6 


1.63 ± 0.74 


1.42 ± 0.91 


0.61 


7 


2.78 ± 1.21 


1.65 ± 1.05 


0.85 


8 


1.82 ± 0.81 


0.94 ± 0.63 


0.61 


9 


2.93 ± 1.39 


3.30 ± 2.11 


0.62 


10 


2.12 ± 0.95 


1.16 ± 0.78 


0.60 


11 


1.04 ± 0.50 


1.94 ± 0.91 


0.48 


12 


7.05 ± 3.25 


3.72 ± 2.69 


0.47 



^ Ljji is the 8-l,000|J.m bolometric infrared luminosity. 

^ The 8-4OM0 SN dust is assumed. For the MW dust, a factor 

of 0.78 should be multiplied. 

"^ r is the radius of dust shell. It can be apprpximated as 

T ~ ^^TxTy which were obtained from HST I8i4 images. 



3 ANALYSIS 

We take a similar approach as performed in our previous 
study. Synthetic spectra of stellar populations reddened by 
various dust models are applied to reproduce the data of the 
12 sample g:ala xies. For the fu ll details of this approach, the 
reader refers to lKawara et alJ (^2011 ) and references therein. 



3.1 Dust 

Our dust models are a combination of an optical properties 
of dust and a geometry of the dust distribution. 



grains distribute both within the system of stars as internal 
dust and in the screen in front of the system. 

The effects of dust also depend on the physical dis - 
tance between the dust and the stars (jCalzetti et al.lll994l V 
If the dust is located physically distant from the stars, the 
effect of dust is to remove photons from the line of sight 
through absorption and scattering. Hereafter, we call this 
distant configuration no scattering. On the other hand, if 
the dust is located close to the stars, the effect is not only 
to remove photons from the line of sight but also to convey 
photons back into the line of sight by scattering. We exam- 
ine this close configuration in two extreme cases, namely, 
forward- only scattering and isotropic scattering. In the in- 
ternal dust geometry where the dust is located close to the 
stars, we account for the effect of the scattering through the 
relations Te// = (1 — 7) Text for forward-only scatte ring and 
Jeff — (1 — 7)°"^^ Text for isotropic scattering (Kaw ara et al.l 
|201l|), where 7 is albedo, Te// the effective absorption op- 
tical depth. Text the sum of the optical depth of absorption 
and scattering. In the foreground screen geometry, if the 
dust screen is physically distant from the stars, no scat- 
tering applies. On the other hand, if the dust screen is 
physically close to the stars, the scattering can put pho- 
tons into the line of sight. We thus consider three cases 
for the foreground geometry: Te// = Text for no scatter- 
ing, Te// = (1 — ^)Text for forward-ouly scattering, and 
^eff = (1 — 7)°"^Tea;t for isotropic scattering. The relation 
between the observed flux f^, and the unobscured flux /° 
is: fvl fv — (X ~ ^~'^''^^)/^eff for an internal dust geometry 
and /^//° = e-^-^ with Tap = A^(l - e'^-ff^^) for a dust 
screen geom etry with the number o f N clumps along the 
line of sight (jNatta k, Panagialll984l V It is noted that there 
is a definite difference between attenuation and extinction 
curves; we regard attenuation curves as fvj jv and extinction 

curves as Text- 



3.1.1 Geometries of dust distribution 

The attenuation produced by dust in front of an extended 
source depends not only on the optical properties of dust 
grains, but also on the geometry of spatial distribution of 
dust, especially when starlight is attenuated by a thick, 
clumpy dust layer or dust is internal of the extended source. 
Analytic solutions to determine the dust p arameters for var- 
i ous g eom etries are ext e nsively stud ied bv lNatta k, Panagial 
(|l984l ) and lCodd (ll973l ). ICodJ(ll973l ) applied the two-stream 
approximation to the radiative transfer for a circumstel- 
lar dust shell and this approximation was developed for a 
dust slab containing emitters and absorbers in Appendix by 
iKawara et al.l (|201lh . It is noted that the slab approach is 
a good approximation to a spherical geometry. The reader 
may also refer to a n excellent review and discussion in 
ICalzetti et~all(|l994l ). 

We consider the following geometries of dust distribu- 
tion: (1) Internal dust geometry - absorbers (dust grains) 
and emitters (stars) are uniformly mixed and distributed; (2) 
Foreground dust screen geometry - absorbers are physically 
separated from and located in a screen in front of emitters, 
where we assume that there are N clumps in average along 
the line of sight, thus the dust distribution is clumpy (small 
N) or uniform (very large A); (3) Composite geometry - dust 



3.1.2 Extinction curves 

We examine two extinction curves synthesized in eje ct a of 
SNe I I, the empirical extinction curve deri ved by iCalzettil 



(|200lh . the MW extinction curve given bv IPraind (|2003 
and the pseudo MW extinction curve. 

The MW extinction curve exhibits a pronounced feature 
at 2175 A. This feature is known to be absent in the typi- 
cal extinction curves of the Small Magellanic Clouds (SMC) 
and some galaxies. To check the effect of the feature on the 
extinction law in our galaxies, we create the pseudo MW 
extinction curve which is made by removing the 2175 A fea- 
ture, namely, the extinction bump and albedo dip, from the 
MW extinction curve. 

The Calzetti extinction curve is based on the attenua- 
tion law that was derived in a range from 0.12 to 2.2 \\jn. by 
analyzing local starbursts and blue compact dwarf galaxies 
(JCalzetti et al.lll99J l . ICalzettil (|200lh expresses this law in a 
functional form combined with a geometry of a uniform dust 
screen with no scattering (i.e., fu = e~^'" f^). This cannot 
be used for longer wavelengths, because the extinction curve 
(tu) sharply drops beyond 2.2 ^tm. The Calzetti extinction 
curve is almost identical to the MW curve from 0.6 to 1.6 
|j,m. We thus extend the Calzetti curve to longer wavelengths 
by adopting the MW extinction curve for wavelengths ^1.6 
\iin. Because the Calzetti curve has already been processed 
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Table 6. Parameters of synthetic template spectra. 



Parameter 



Range 



Z 0.005 - 2.5Z0 at 0.005,0.02,0.2,0.4,1.0,2.5^0 
t 1 - 6,300 Myr in a step of 0.2 dex 

tsfr 10 - 10,000 Myr in a step of 0.3 dex 
^o.3|J.m 0-15 mag in a step 0.1 mag 

15 - 100 mag in a step of 1 mag 



Note. - Z, t, tsfr-) ctnd Ao.3|j.in are the metaUicity, age after the 
onset of the star formation in Myr, e- folding timescale of the 
star formation in Myr, and extinction at 0.3|J.m in magnitude, 
respectively. 



spectra are generated in the parameter range as shown in 
Table [1 

The equivalent widths of [OII]A3727 and MgIIA2798, 
and D(Call) were calculated from the synthetic spectra. 
[Oil] is a nebular line and unpredictable directly from the 
synthetic spectra. We first obtain the Hq, luminosity by 
counting ionizing photons in the synthetic spectra. The 
metallicity scaling of [Oil] /Ha f or the ionized param eter 
^ = 3 X lO'^cm s~^ was given bv iKewlev et al.1 (|2004l V We 
thus adopt their theoretical curves with an additional con- 
dition, which is applied to avoid a negative or unrealistic 
value, of [OII]/Ha > O.l. lKewlev et al.1 (|200 J ) estimated the 
uncertainty in the [Oil] intensity to be ^-^30%. 



by the foreground screen, we consider this geometry only for 
the Calzetti curve. 

Combining the optical constants in the literature with 
the grain s pecies and size distribution of SN-condensed dust 
derived by Nozawa et al.l ([2003. ) . and using the Mie theory, 
iHirashita et al.l (|2005|) calculated absorption and scattering 
cross sections of homogeneous spherical grains with various 
sizes condensed in mixed and unmixed eject a. Comparing 
their models to the exti nction curve of SPSS J 1048+637 (a 
BALQSO at z - 6.2), IHirashita eraD (|2005l ) suggest the 
mass range of SN progenitors to be < 30 Mq, whereas the 
upper mass of SN progenitors is usually assumed to be 40 
y^c ^- We thus adop t the two extinction curves of SNe II for 
the ISalpeterl (|l955l l Initial Mass Function (IMF) ; one has a 
mass range of progenitors from 8 to 30 Mq , and the other 
from 8 to 40 Mq. These were obtained by interpolating 
or ex trapolating the extinction curves that IHirashita et al.l 
(|2005f) calculated for unmixed SNe II with progenitor masses 
of 13, 20, 25, & 30 Mq. We adopt these in the present 
study and call the original SN extinction curves in a sense 
that the subsequent processing of dust destruction are not 
considered. Recent studies indicate that the original SN ex- 
tinction c urves become flatter through the reverse shock de - 
struct ion ([Bianchi k SchneiderllJOOTh IHirashita eraIll20Q8h , 
while IHirashita et al.l (|2010l ) shows the importance of the 
shattering in a warm ionized medium (WIM) in starburst 
environments and this effect makes the extinction curves 
steeper. They suggest that the final extinction curves, af- 
ter the reverse shock destruction and the shattering, have 
a slope similar to that of the original curve. The degree of 
the reverse shock destruction and shattering critically de- 
pends on the assumed parameters, especially gas density of 
the ISM. Considering these uncertainties and the similarity 
between the original and final extinction curves, we decided 
to adopt the original SN extinction curves in the present 
study. 



3.2 Star formation history 

The evolutionary synthesis codes by iBruzual k, CharlotI 
(|2003|) are used to generate synthetic spectra of evolution- 
ary stellar population models. We adopt the Salpeter (1955) 
IMF with the mass cutoffs of O.l-IOOM©. The star forma- 
tion histories are exponentially declining models with the 
star formation rate (SFR) oc t~i exp(— t/ts/r), where tsfr 
is the e-folding timescale and t is the age. The synthetic 



3.3 Minimum x^ method 

We search for the best set of parameters by using a minimum 
X^ method as defined below: 



X 



■ XSED + ^(X[OII] + XEW(abs) + XLj^ + Xm^) 



il) 



The weighting factor W is introduced to balance the 
weight between the SED datsfj- 
Each term is defined as: 



2 
XSED 



^^ filters 



- hFtemp) /CTSED^ 



(2) 



X[OII] = {EW{[Oll])o,s - EW{[0ll])tempf/c7loil]. (3) 
XEW(abs) = {EW{Bhs)ohs - EW{Bhs)tempf/(JEW{^hs) (4) 



2 
XLjji 

XMd 



{LiR — LiR temp) /cTLjf 
■ Md temp) /^M^' 



(Md 



(5) 
(6) 



Equation [2] is related to photometric data of the contin- 
uum from stellar atmosphere. 6 is a normalization constant 
between the template fiux density Ftemp and the observed 
Fohs^ and the number of photometric filters Nfuters from UV 
to 5.8 \im. 8.0^im data are excluded from the SED data, be- 
cause dust emission may be significant in the 8.0^im band. 
In order to avoid too small and non-realistic photometric 
error, 5% of the fiux is added in the quadratic form to the 
photometric error aph, i.e., (Tsed = o-ph + (O.ObFobs)'^ • 

Equation [3] is related to the [Oil] line from HII re- 
gions ionized by you ng OB stars. EW{[OII]) predicted using 
iKewlev et al.l (200J) theoretical curve contains ^30% uncer- 
tainty, thus a^ou] = ^%waou])obs + (0-3^V^([OII]),,^p)^ 
where <^EW([oii])obs i^ the observed error. Equation |4] mea- 
sures the absorption lines in stellar atmosphere. EW{abs) is 
EW{MgII) for HRIO and D{CaII) for the others. cr|vv.(abs) 
is the observed error. 

Equation [5] compares Ljr given in Table [5] with the 
prediction from the model, i.e., Ljr temp oc /(/^ — fu)dh', 
where /° and fu are intrinsic and reddended synthetic spec- 
tra. Equation [6] compares Md derived from fir-IR data with 



Md 



temp 



the dust mass estimated from the dust extinction 



^ The number of photometric data points was chosen so that 
the SED fitting did not overweight the other data. However, we 
found at a late stage of the analysis that the weight of the SED 
still overwhelmed the other data for many galaxies. W was thus 
introduced to ensure that the x^ of the SED has the similar mag- 
nitude as that of the other data. As a result, we obtained W = 1 
for galaxies 2, 4, 6, 12, and W = 3 for the others 
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(i.e., optical depth) and the assumed dust geometry. Md temp 
was derived as follows. We first assume a composite geom- 
etry consisting of screen dust and internal dust. Suppose a 
spherical shell with a radius of r and a thickness of 5r^ where 
the internal and screen dust fills inside and in the shell, re- 
spectively. Then, assuming 5r/r <^ 1, the mass of the inter- 
nal dust is Min = 47r(rin/^ext)^^/3, while that of the screen 
dust is Msc = ^7r{rsc/i^ext)r'^- nn and Tsc are Text for inter- 
nal dust and screen dust, respectively, and they are related 
to r and 6r as nn = Pini^extr and Tsc = psd^extSr, where pin 
and psc are the mass density of the internal and screen dust, 
respectively. The dust mass is derived directly from nn, Tsc, 
and r with no needs to assume 6r, pin and psc- r can be ap- 
proximated as r ^ yJ^xTy which were obtained from HST 
/814 images, where Vx and Vy are the semi-major and the 
semi-minor axes, respectively, r is given in Table [5] 



4 RESULTS 

The results of the minimum x^ analysis are illustrated in 
Figure m The left panels of the first 3 rows plot x^ (=X^/^5 
where h> is degree of freedom jj as a function of nn/ijin+Tsc) 
for a composite geometry consisting of uniform screen dust 
and internal dust; nnjijin + Tsc) = for the foreground 
screen geometry, while nn/{nn + Tsc) — 1 for the internal 
dust geometry. The scattering properties are assumed to be 
no scattering, isotropic scattering (^ = 0) , or forward-only 
scattering (^ = 1). In case of no scattering, the scattering 
properties are no scattering for screen dust and forward scat- 
tering for internal dust because no scattering cannot be as- 
sumed for internal dust. The middle and right panels of the 
first 3 rows show xt as a function of N the number of clumps 
along the line of sight in the foreground screen geometry and 
the composite geometry with nn/{nn -\-rsc) = 0.4. The sev- 
eral linestyles are used to specify the extinction curve: Solid 
lines to the 8-30 Mq SN curve, dash to the 8-40 Mq SN, 
dot to the Calzetti, dash dot to the MW, and dash dot dot 
dot to the pseudo MW. 

The four panels in the forth row plot Ao.s^im the ex- 
tinction (nn + Tsc) in magnitude at 0.3^Lm, t the age (i.e., 
the time after the onset of the star formation) in Myr, Z the 
metallicity in solar units, and tsfr the e-folding timescale of 
star formation in Myr as a function of nn/inn+Tsc)- The left 
panel in the bottom row plots the best-fit SED along with 
the unreddened SED and the right panel the contour map 
of xl ~ X?(^^^) where 'x^{'min) is the smallest xi among 
all the models. 

The trends of the results are:(l) the SN dust curves gen- 
erally have reducing x^ as nn/(Tin + Tsc) is approaching to 
zero (more dominated by screen dust) and N approaching to 
oo (more uniform screen); (2) Some galaxies have a small xi 
in a clumpy foreground screen geometry, especially, galaxies 
Nos. 5 and 12 have the minimum xS at N of 10 and 3, respec- 
tively; (3) The MW and pseudo MW curves sometimes have 
reducing x^ as a geometry gets more dominated by internal 



"^ Degree of freedom is defined as ?/ = ^/^^ter + 4VK — 5 where 
^filter is the number of photometric data used for Xsed ^^^ ^ 
the weighting factor. The number of fitted parameters is 5 (i.e., 
b, t, tsfr, -4o.5|u.m, Z); Nfuter = 9 for the COSMOS galaxies, 13 
for SST J1604, and 10 for HR 10. 



dust at nn/{nn + Tsc) = 0.8 - 0.9. The data point around 
rest 4|j,m in the SED plots is the IRAC 8.0 ^im data, and 
shows an excess from the model for most of the galaxies. As 
discussed in section 3.3, the data of this band were not used 
in the analysis, because dust emission may be significant in 
this band. The SED plots in Figure [H indeed suggest dust 
emission around 4|j,m. 

Table [71 compares the smallest xl of the respective ex- 
tinction curves along with the dust geometry, the scatter- 
ing property, and the total extinction including internal and 
screen dust. There are two obvious features: (1) the opti- 
cal depth (Ao.3p,m) in the geometry dominated by internal 
dust is several times greater than that in the foreground 
screen geometry; (2) the MW curve generally has the small- 
est xt in the composite geometry dominated by internal 
dust, whereas the SN curves do in the foreground screen 
geometry, and the pseudo MW dust model is in midway. 
Finding of an higher optical depth for the internal absorp- 
tion is an expected and obvious result, because the screen 
geometry always maximizes the extinction. 

The contour plots in Figure [H appear to show the 
age (t) /extinction (Ao. Sum) degeneracy, suggesting that t and 
tsfr given in Tableware not very accurate. The uncertain- 
ties in the age determination were propagated to the uncer- 
tainty to determine the stellar mass of galaxies, because the 
luminosity to mass ratio strongly depends on the age of stel- 
lar populations. This degeneracy can be resolved by adding 
data from lOO^im to the submillimeter, greatly reducing the 
present uncertainty in Ljr. 

Table [8] summarizes the best-fit models. Out of our 12 
galaxies, 10 (83%) have a pure foreground dust screen, and 
9 of the 10 galaxies have a uniform foreground screen. 8 
(67%) of the 12 galaxies have the best-fit model with one of 
the SN dust curves, all of which are associated with a pure 
foreground dust screen. Internal dust is only significant in 
two galaxies, all of which have the best-fit models with the 
MW or pseudo MW extinction curve. Because the extinc- 
tion curves analyzed here can be ordered from the steepest 
to the flattest as the MW, pseudo MW, Calzetti, 3-40 Mq 
SN, and 3-30 Mq SN, the MW and pseudo MW curves are 
too steep to agree with a pure foreground screen, and need to 
be associated with a geometry dominated by internal dust or 
a clumpy foreground screen to make the attenuation curve 
flat or gray. This requires a few times more dust mass than 
that in a pure uniform foreground screen. Our results indi- 
cate that the 8-40 Mq SN curve with a uniform foreground 
screen is the best dust model to account for z ^ 1 ULIRGs. 



5 DISCUSSION 



I Gordon et al.l (|l997l l applied models of stars and dust to the 
SEDs of 30 local starburst galaxies using the MW and SMC 
extinction curves, and found that the dust has an extinc- 
tion curve lacking a 2175 A feature with a slope interme- 
diate between the MW and SMC curves. Th e dust geome- 
try has a pure foreground screen with clumps. iGaskell et al.l 
(;2004) found indications of a very flat extinction curve in 
the ultra-violet (UV) from an analysis of 1,090 AGNs. Large 
samples of SDSS quasars indicate that the extinction law in 
quasars at z < 2.2 is described by S MC reddening but not 
by reddening in the LMC and MW JRichards et all l2003l : 
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Figure 4. Results of the minimun x^ analysis. The source number and the scattering property are indicated in each panel. The left 
panels in first 3 row show x? for a uniform screen as a function of Tin/{Tin + Tsc)- The middle and right panels in first 3 row represent 
a screen of A^ clumps for Tinlijin + Tsc) = and Tin/(Tin + Tsc) = 0.4, respectively. In case of no scattering, the scattering properties 
are no scattering for screen dust and forward scattering for internal dust. Solid lines represent 8-30 Mq SN dust, dash 8-40 Mq SN 
dust, dash dot MW dust, and dash dot dot dot pseudo MW dust, dot the Calzetti dust. The ^t/i row panels show Tin + 'T'sc at 0.3|J,m in 
magnitude Ao.3|j,m, the metallicity Z, the age t(yr), and the timescale tgfr as a function of Tin/{Tin -\-Tsc)- The bottom row panels show 
the best-fit SED along with the unreddended SED and the contour map of xt ~ X?(^^^^)- The full version of the paper is for online 
publication only, where the full version of this figure is available. 
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Table 7. Comparison of dust models 



No. 




Reduced 


X square (x?") 






Scattering property 








Tin/ 


Tin-\-sc 


{N^) 






Ao.SM-m 








MW 


pMW 


Cal 


SN30 


SN40 


MW 


pMW 


Cal 


SN30 


SN40 


1 


2.18 


2.72 


2.18 


5.27 


0.79* 


Isotrop 


Isotrop 


- 


No 


No 




0.9 (-) 


0.0 (oo) 


- 


0.0 (oo) 


0.0 (oo 


) 17.5 


7.1 


4.3 


2.4 


4.3 


2 


3.35 


2.87 


2.04 


3.65 


1.61* 


Isotrop 


No 


- 


Isotrop 


No 




0.9 (-) 


0.6 (-) 


- 


0.0 (oo) 


0.0 (oo 


) 14.4 


7.3 


3.9 


4.1 


3.4 


3 


2.09 


1.98 


2.42 


2.51 


1.62* 


Isotrop 


Isotrop 


- 


No 


No 




0.8 (-) 


0.0 (oo) 


- 


0.9 (-) 


0.0 (oo 


) 23.8 


8.8 


5.9 


18.8 


6.0 


4 


3.93 


0.86* 


3.02 


8.72 


2.33 


Forward 


Isotrop 


_ 


No 


No 




0.9 (-) 


0.0 (oo) 


- 


0.0 (oo) 


0.0 (oo 


) 25.0 


5.9 


2.8 


2.5 


3.5 


5 


1.52 


2.19 


2.07 


2.45 


0.80* 


No 


No 


_ 


No 


No 




0.9 (-) 


0.7 (-) 


- 


0.0 (oo) 


0.0 (10 


) 13.1 


7.5 


3.3 


2.9 


4.8 


6 


1.88 


1.72 


1.61 


1.51 


1.37* 


No 


Isotrop 


- 


No 


No 




0.8 (-) 


0.0 (oo) 


- 


0.0 (oo) 


0.0 (oo 


) 11.9 


6.8 


4.3 


2.8 


4.9 


7 


4.30 


4.59 


5.00 


5.61 


2.35* 


No 


Isotrop 


_ 


No 


No 




0.0 (10) 


0.0 (oo) 


- 


0.0 (oo) 


0.0 (oo 


) 3.3 


4.4 


3.0 


2.6 


3.6 


8 


1.35 


1.99 


1.51 


3.00 


0.98* 


No 


Isotrop 


- 


Isotrop 


No 




0.0 (10) 


0.0 (oo) 


- 


0.0 (oo) 


0.0 (oo 


) 4.0 


5.4 


3.6 


2.9 


3.6 


9 


2.15 


1.02* 


1.51 


2.10 


1.88 


No 


Isotrop 


- 


No 


No 




1.0 (-) 


0.9 (-) 


- 


0.0 (oo) 


0.0 (10 


) 32.5 


21.3 


3.5 


3.6 


4.9 


10 


2.88 


1.23 


0.75* 


2.37 


2.08 


Isotrop 


Isotrop 


- 


Foward 


No 




0.8 (-) 


0.0 (oo) 


- 


0.0 (oo) 


0.0 (oo 


) 9.4 


4.5 


3.6 


4.0 


3.6 


11 


5.09 


5.55 


5.05 


3.56* 


3.75 


No 


Isotrop 


_ 


No 


Foward 




0.0 (10) 


0.9 (-) 


- 


0.0 (oo) 


0.0 (oo 


) 8.8 


42.5 


5.8 


5.4 


17.5 


12 


0.84* 


0.84* 


1.36 


1.71 


1.41 


Isotrop 


Isotrop 


- 


No 


No 




0.4 (3) 


0.4 (3) 


- 


0.0 (oo) 


0.0 (oo 


) 20.0 


20.0 


5.0 


4.0 


4.4 



Note. - The smallest xJ is given per galaxy per dust model. "*" marks xJ of the best-fit model for each galaxy. 'MW, 'pMW, 

'Car, 'SN30', and 'SN40' denote a dust model of the MW, pseudo MW, Calzetti, 8-3OM0 SN, and 8-4OM0 SN, respectively 

'Isotrop', 'Forward', and 'No' denote isotropic, forward, and no scattering, respectively. 

" X? = X?/^5 where u is degree of freedom. 

^ The number of clumps N is given in parentheses for a geometry dominated by foreground screen dust. 



iHopkins et HI l2004l ) . iGallerani 6^11 (|2010l ^ proposed that 
reddening of quasars at z — 3.9 — 6. 4 has a tendency to be 
fiat at A < 2000 A. Mai olino et~all ([2004) found that the 
BALQSO at ^=6.2, requires a flat extinction curve which 
is quite different from those observed in the MW, LMC 
and SMC, but in excellent agree ment with the synthetic 
SN curves (|Todini fc Ferraral[200lh . The flat extinction law 
that we found m. z ^ 1 ULIRGs are probably incompatible 
with the steep SMC-like extinction law reported for z <2.2 
AGNs, and in agreement with the flat extinction law found 
in high redshift quasars. 

It should be reminded that the SN-dust models gener- 
ally predict flat extinction c urves but there a re int rinsic dif- 
ferences among the mo dels. Maiolino et al.l (120041 ) adopted 



the SN-dust model by iTodini k Ferraral (l200lh to derive 
the extinction curves expected in high-z quasars. The cal- 



culated curves show a flat feature at UV wavelengths as a 
result of the main contribution of amorphous carbon grains 
with radii around 300 A, added by the minor contribu- 
tions of small silicate and FesO^ grains with radii around 
20 A. lBianchi fc Schneide^ (l20(37l) rev isited the dust forma- 
tion model bv lTodini k, Ferraral (|200l|) and also followed the 
evolution of the dust grains from the time of formation to 
their survival through the passage of the reverse shock. They 
showed that the reverse shock predominantly destroys small 
silicate and Fe3 04 grains, and the resultant extinction curves 
become flatter, being dominated by carbon grains whose ex- 
tinction coefficients are rela tively insensitive to wavelengths. 
On the other hand, jHirashita et all (|2005h calculated the 
extinc tion curves using SN Il-dust model bv iNozawa et al.l 
(|2003l V The derived extinction properties more or less de- 
pend on the SN progenitor mass, but are quite flat as a re- 
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Table 8. Properties of the best-fit models. 



No. 


Tin/Tin+sc (N) 


Dust model 


Scattering 


xl 


Z 


t 


tsfr 


-4o.3M.m 


M,^ 


Msc 


M* 


Md/SN 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


1 


0.0 (oo) 


SN40 


No 


0.79 


0.2 


1500 


10000 


4.3 


_ 


1.6 


27.5 


0.0079 


2 


0.0 (oo) 


SN40 


No 


1.61 


1.0 


800 


630 


3.4 


- 


1.9 


6.5 


0.031 


3 


0.0 (oo) 


SN40 


No 


1.62 


1.0 


60 


10 


6.0 


- 


1.0 


1.4 


0.16 


4 


0.0 (oo) 


pMW 


Isotrop 


0.86 


1.0 


100 


40 


5.9 


- 


1.3 


1.2 


0.19 


5 


0.0 (10) 


SN40 


No 


0.80 


0.4 


60 


20 


4.8 


- 


2.8 


1.9 


0.25 


6 


0.0 (oo) 


SN40 


Isotrop 


1.37 


1.0 


1500 


320 


4.9 


- 


1.7 


11.3 


0.017 


7 


0.0 (oo) 


SN40 


No 


2.35 


0.4 


40 


10 


3.6 


- 


2.5 


0.6 


0.38 


8 


0.0 (oo) 


SN40 


No 


0.98 


1.0 


150 


40 


3.6 


- 


1.2 


1.5 


0.085 


9 


0.9 (-) 


pMW 


Isotrop 


1.02 


1.0 


150 


160 


21.3 


2.4 


0.8 


1.5 


0.23 


10 


0.0 (oo) 


Calzetti 


No 


0.75 


0.4 


60 


10 


3.6 


- 


1.0 


1.0 


0.12 


11 


0.0 (oo) 


SN30 


No 


3.56 


1.0 


630 


160 


5.4 


- 


0.9 


3.6 


0.058 


12 


0.4 (3) 


MW 


Isotrop 


0.84 


0.4 


2000 


630 


20.0 


0.7 


2.2 


42.0 


0.0098 




0.4 (3) 


pMW 


Isotrop 


0.84 


0.4 


2000 


630 


20.0 


0.7 


2.2 


42.0 


0.0098 



Note. - Columns 1-5 are defined in Table O 

Columns 6-9: Metallicity(Z) in solar, age(t) in Myr, e-folding timescale(ts/^) in Myr, and extinction at 0.3|xm(Ao.3|xm) in 

magnitude, respectively. 

Columns 10-11: Mass of internal dust (Min) and screen dust (Msc), respectively, in units 10^ solor mass. 

Column 12: Stellar mass (M*) in units 10-*^ -"^Mq that is estimated as Ljji/(L/M^), where L/M^ is the luminosity to mass ratio for 

the stellar population models. 

Column 13: Dust mass(Mc/) per SN where (M^) in solar units. The number of SNs per galaxy is derived from the IMF of the 

model of star formation history and integration of all stars with > 8 Mq. 



suit of the dominant contribution of large (> 0.1 ^im) amor- 
phous carbon or Si grains. As discussed in Se ction 3.1.2, the 
rever se shock efficiently destroys small grains (jNozawa et al.l 
[20071), making the extinction curves mu ch flatter than the 
original curves fsee lHirashita et al.ll2008l ). Then large grains 
that survive the destruction through the reverse shock are 
shattered in the WIM, producing a large number of small 
grains. This finally leads to the ex tinction curves simila r 
to the original SN extinction curves (jHirashita et al.ll201Ql ). 
The growth of dust grains through acc retion in the ISM 
would also be efficient in infra red galaxies (|Dweklll998llDrainel 
l2009llMichaIowski et al.ll20ld ) , and also cause flat extinction 
curves. 

As shown in Table [51 our z ^ 1 ULIRGs have an age 
of mostly t ^ 650 Myr with four exceptions (i.e., t = 800, 
1,500, 1,500, 2,000 Myr), a stellar metalli city of ^ = 0.2-1.0 
solar, and a dust mass of Md ^ lQ^M0.'MichaIowsk i et al.[ 
(|201Q| ) suggested that AGB stars are not sufficient enough 
to form dust at such an early epoch of galaxy evolution his- 
tory, whereas SN-condensed dust followed by grain growth 
in the ISM is important to account for the inferred dust 
mass. It is noted that our analysis cannot rule out a possi- 
ble contribution from AGB stars, becau se some observations 
(|Groenewegen[[l997[| Ganger et al.[[ 19991 ) suggest that large- 
sized dust grains are formed in AGB stars and larger grains 
have flat extinction curves. The dust mass per SN in our 
z ^ 1 ULIRGs mostly ranges from 0.01 Mq/SN to 0.3 
Mq/SN. This should be compared with theoretical predic- 
tion of the dust yield per SN, which is 0.01 - 0.1 Mq/SN 
with dust destruction and 0.1-1 Mq/SN without dust de- 
struction. If the ISM process of dust growth is important, 
the expected yield will be more than these. However there 
would be a possible bias effect on the dust mass per SN 
inferred from the observations, because a significant corre- 
lation between the age t and Ma/ SN can be seen in Table 



[8] This correlation seems spurious and simply reflects the 
uncertainty of determining t caused by the t and Ao.sp^m 
degeneracy, because the mass-to-luminosity ratio strongly 
depends on t. 

Our results indicate that the dust geometry of a pure 
foreground dust screen generally explains the data of the 
z ^ 1 ULIRGs better than that dominated by internal dust. 
Color-color plots of optical and UV data of local starbursts 
were compared with various dust distribution models by 
[Gordon et al.[ (|l997[), who found that the observed spread 
cannot be explained by the internal dust geo metry, but the 
foreground screen geometry can explain it. [ Calzetti et al.[ 
(1994) also showed the failure of the internal dust geome- 
try to explain the spread of data of local starbursts, sug- 
gesting the foreground screen geometry with a fiat extinc- 
tion curve (i.e., Calzetti curve). [Meurer et al.[ (|l995[ ) showed 
that a strong correlation between the far-IR excess and the 
UV spectral slope for local starbursts can be explained by 
dust in a foreground dust screen geometry with the Calzetti 
curve, but not by an internal dust geometry. They attributed 
the dust in a foreground screen to the Galactic winds fre- 
quently observed in starbursts that will sweep out any dif- 
fuse ISM from the starbu rst sites on a timescale of a few 
Myr (JHeckman et al.[[l99Q[ V The foreground dust screen ge- 
ometry we are looking at in 2; ^ 1 ULIRGs may also result 
from the dust that was swept by the Galactic winds. 



6 CONCLUSIONS 

In order to study the extinction law in galaxies (ULIRGs) at 
2; ^ 1, we have analyzed the multi- wavelength photometric 
and spectroscopic data of 12 ULIRGs which are archived or 
published in the literature; 10 galaxies come from the COS- 
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MOS survey catalogs, and the others are SST J1604+4304 
and HRIO. 

The data are compared with models of stars and dust. 
We applied five extinction curves, namely, the Milky Way 
(MW), the pseudo MW which is MW-like without the 2175 
A feature, the Calzetti, and two SN dust curves, by combin- 
ing with various dust distributions, namely, the foreground 
screen geometry, the internal dust geometry, and the com- 
posite geometry with a combination of dust screen and in- 
ternal dust. Several scattering properties are used: no scat- 
tering, isotropic scattering, and forward scattering. In case 
of no scattering the effect of dust to convey photons into 
the line of sight is ignored, while this effect is included in 
isotropic scattering and forward scattering. 

The two SN extinction curves tha t we applied here 
were calculated bvi Hirashita et al.l (|20Q5l ) bas ed on the g:rain 
Tope rties of SN-condensed dust derived by iNozawa et al.l 
2003r ). These SN extinction curves were calculated assum- 
ing the unmixed ejecta of SNe II for the Salpeter IMF in 
a progenitor mass range of 8 - 30 or 8 - 40 Mq and do 
not include the effects of the reverse shock in the ejecta and 
interaction between SN dust grains and the ambient ISM. 

Employing a minimum x^ method, we find that the 
foreground dust screen geometry, especially combined with 
the 8-40 Mq curve, provides a good approximation to 
the real dust geometry, whereas internal dust is significant 
only for 2 galaxies. The SN extinction curves, which are 
flatter than the other extinction curves, reproduce the data 
of 8(67%) galaxies better than the other curves do. Our z ^ 
1 ULIRGs require an age of mostly t ^ 650 Myr, a stellar 
metallicity of z = 0.2 — 1.0 solar, and a dust mass of Md ^ 
10^ Mq. Considering these physical features, we conclude 
SN-origin dust is the most plausible to account for the vast 
amount of dust masses and the flat slope of the observed 
extinction law. The inferred dust mass per SN ranges from 
0.01 to 0.4 Mq/SN. 



7 ONLINE MATERIALS 

The full version of the paper is for online publication only, 
where the full version of Figure U is available. 
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